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ABSTRACT 

We report the first results of a programme aimed at studying the properties of 
high redshift galaxies with on-going massive and dominant episodes of star formation 
(HII galaxies). We use the i(H/3) — a distance estimator based on the correlation 
between the ionized gas velocity dispersions and Balmer emission line luminosities 
of HII galaxies and Giant HII regions to trace the expansion of the Universe up to 
z ^ 2.33. This approach provides an independent constraint on the equation of state 
of dark energy and its possible evolution with look-back time. 

Here we present high-dispersion (8,000 to 10,000 resolution) spectroscopy of HII 
galaxies at redshifts between 0.6 and 2.33, obtained at the VLT using XShooter. 
Using six of these HII galaxies we obtain broad constraints on the plane Q,m — wq. 
The addition of 19 high-z HII galaxies from the literature improves the constraints 
and highlights the need for high quality emission line profiles, fluxes and reddening 
corrections. The 25 high-z HII galaxies plus our local compilation of 107 Hii galaxies 
up to z = 0.16 were used to impose further constraints. Our results are consistent 
with recent studies, although weaker due to the as yet small sample and low quality 
of the literature data of high-z HII galaxies. 

We show that much better and competitive constraints can be obtained using a 
larger sample of high redshift HII galaxies with high quality data that can be easily 
obtained with present facilities like KMOS at the VLT. 
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1 INTRODUCTION 

H II galaxies (HUG) are most extreme narrow emission line 
star forming systems selected from spectroscopic surveys as 
those having the largest equivalent widths (EW) in their 
emission lines, i.e. EW(H/3) > 50A (or EW(Ha) > 200A) 
in their rest frame and being extremely compact. The lower 
limit required in the selected equivalent width of the re- 


* Partly based on observations obtained at the European South¬ 
ern Observatory, programme ID 091.A-0413(A) 
f E-mail: rjt@inaoep.mx & rjt@ast.cam.ac.uk 


combination hydrogen lines is of fundamental importance to 
guarantee that the sample of HIIG is composed by systems 
in which a single and very young starburst, less than 5 Myr 
in age, dominates the total luminosity output. This selection 
criterion also minimizes the possible contamination by an 
underlying older population or older clusters inside the spec¬ 
trograph ape rture (cf. [Melnick, Terlevich fc Terlevich|2000 
|Dottori]1981| [Dottori fc Bica|1981| [Chavez et al.|2014[ ) and 

the ionizing photon escape. HIIG thus selected are spectro¬ 
scopically indistinguishable from young Giant Extragalactic 
Hii Regions (GEHRs) in nearby galaxies (e.g. 30 Doradus 
in the LMC or NGC 604 in M33). This is highlighted by the 
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fact that when discovered, the prototypical HIIG I Zwl8 and 
II Zw40 were named “Isolated Extragalactic Hll regions” 
( Sargent fc Searle]|1970| . On the other hand, the compact¬ 
ness requirement biases the sample towards single bursts 
with relatively small crossing times. 

Under the working hypothesis of instantaneous coeval 
star formation (ISF) it is possible to define an age for the 
stellar population. Even if continues star formation (CSF) 
is assumed and a wide range of stellar ages is present the 
ionization will be dominated by the stars born in the last 
5 Myrs. Thus from this point of view not much difference 
is expected in the observables of these extreme emission 
line starbursts between the instantaneous and the contin- 
uos case. There are parameters however that will differ like 
the EW of the IR-Call triplet or the total L(Ha/Mass re¬ 
lation. Another important difference is that in the case of 
CSF one expects to detect the WR features in all systems 
while for ISF they are expected only in a fraction of them. 

The optical properties of HIIG can be considered as 
those of a “naked” extremely young and compact burst of 
star formation rather than those of the host galaxy. 

It has been shown that GEHR and HIIG exhibit a tight 
correlation between the luminosity and the width of their 


emission lines, the L(H/3) — a relation (Terlevich & Melnick 


1981). The scatter in this relation is small enough that it 


can be used to determine cosmic distances independently of 


redshift (see 

Terlevich & Melnick 1 1981 |Melnick, Terlevich & 

Moles 19881 

Melnick, Terlevich & Terlevich||2000| |Fuentes- 


Masip^et^^ljboOO Telles|200^ Bosch, Terlevich & Terlevich 


2002 Siegel et al. 2005| Bordalo &i Telles|201H Chavez et al. 


2012 2014). HIIG can reach H/3 Inminosities larger than 
10“*“ erg s“^ making them observable even at relatively large 
redshifts {z > 3) with present-day NIR spectrographs. 

In (CSF) we showed that the Z/(H/3) — a relation consti¬ 
tutes a viable alternative to SNe la for the determination of 
cosmological parameters and presented a general strategy to 
use high z HIIG as effective cosmological probes to reduce 
significantly the parameter space of the dark energy equa¬ 
tion of state and test its possible evolution with redshift. 

To date, the cosmic acceleration has been traced di¬ 
rectly only by means of SNe la and up to redshift 2 ~ 1.5 


(|Riess et al.||1998[ 

Perlmutter et al.||1999| |Amanullah et al.| 

2010 Hicken et al.j 

2009 Suzuki et al. 12012 ), a fact which im- 


plies that it is important to use alternative probes at higher 
redshifts in order to verify the SNe la results and to ob¬ 
tain more stringent constraints in the cosmological parame¬ 
ters solution space, with the ultimate aim of discriminating 
among the various theoretical alternatives which attempt to 


explain the accelerated expansion of the Universe (cf. Suyu 
et al.|2012| ). 


In this paper we demonstrate the feasibility of using 
HIIG as competitive cosmological tracers to high -2 thanks 
to the availability in 8m class telescopes of highly efficient 
near infrared high spectral resolution spectrographs. 

We would like to reinforce what we consider a central 
point. The fact that a L(H/3) —ct relation exists at all, means 
that it can be used, empirically, as a distance estimator up to 
large redshifts, regardless of the physics that causes it. This 
aspect has been discussed in the literature (among others 


Terlevich & Melnick|1981||Melnick, Terlevich k, Moles 

1988 

Tenorio-Tagle, Munoz-Tunon & Cox 19931 Terlevich 

1997 

Melnick, Tenorio-Tagle & Terlevich|19991 Melnick, Ter 

evich 


fc Terlevich|2000 Zaragoza-Cardiel et al.|2015 ) and will be 


discussed still further as more and better data becomes avail¬ 
able. 


In §2 we present the observations and data reduction, 
and the data (both ours and from the literature) are anal¬ 
ysed in §3. The results are discussed in §4 and conclusions 
and plans for future work are given in §5. 


2 OBSERVATIONS AND DATA REDUCTION 


High spectral resolution spectroscopic observations were ob¬ 
tained using the XShooter spectrograph (Vernet et al. 
at the Cassegrain focus of the ESO-VLT (European South¬ 
ern Observatory Very Large Telescope) in Paranal, Chile 
during the nights of 29 and 30 September 20L|^ We ob¬ 
tained spectra for the three arms using a 0.6 arcsec slit, 
the typical spectral resolution on the VIS arm was ~ 10000, 
whereas in the NIR arm it was ~ 8000. Total exposure times 
per object ranged between 1 and 3 hours. 

A sample of 9 star forming galaxies with rest frame 
equivalent width ( EW) of Yip > 50 A or EW(Ha ) > 200 A 
was selected from |Hoyos et ^ ( |2005[ ) ; |Erb et ( |2006a|b] ) 
and [Matsuda et al. ( |2011| ). The redshift range covered was 
0.64 ^ 2 ^ 2.33. The observation parameters are detailed in 


Table [T] 

A narrow slit was needed to achieve the required reso¬ 
lution but when combined with the excellent seeing resulted 
in pointing problems associated with the fact that these are 
faint and almost stellar objects. All pointings were blind and 
we had to rely on our differential astrometry. This proved to 
be correct, i.e. inside 0.2 arcsec in all cases but one. The nar¬ 
row slit has another negative effect in the sense that slight 
offsets or seeing becoming worse than the slit size result in 
light losses. These aspects highlight the importance of spec¬ 
troscopy with integral field units that minimizes the astrom¬ 
etry problems while at the same time optimizing the light 
collection. 

The data reduction was carried out using the XShooter 
pipeline V2.3.0 over the GASGANO V2.4.3 environmenl[^ 
using the ‘physical model mode’ reduction. 


3 DATA ANALYSIS 
3.1 Emission line widths 

To determine the FWHM of the emission lines, single 
gaussians were fitted to the ID spectral profiles of the 
[O HI] A5007A and Ha lines when availably These fits were 
performed using the IDL routine gaussf it. Figure [37T] shows 
the fits to the [O III] A5007A line. 

The uncertainties of the measured FWHM were esti¬ 
mated using a Montecarlo analysis. A set of random realiza¬ 
tions of every spectrum was generated using the data pois- 
sonian Icr 1-pixel uncertainty. Gaussian fittings for every 
synthetic spectrum in the set were performed afterwards, 

^ Observing programme ID 60.A-9022(C) 

^ GASGANO is a JAVA based Data File Organizer developed 
and maintained by ESO. 

® Historically, H/3 has been used but, when available, the stronger 
Ha line is preferred 
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Figure 1. Gaussian fits to the [O III] A5007A line for the nine observed objects as labelled. Upper panel: The single gaussian fit is shown 
with a dashed line (thick black). The parameters of each fit are shown in the top left corner. Lower panel: Residuals from the fit. 


and we obtained a distribution of “synthetic” FWHM from 
which the Icr uncertainty for the widths measured in the 
spectra follows. 

Table lists the FWHM measurements for the high 
resolution observations prior to any correction such as in¬ 
strumental or thermal broadening. Column (1) is the target 
name, column (2) is the heliocentric redshift as measured 
from the observed emission lines, columns (3) to (6) con¬ 
tain the measured central wavelength in nm and FWHM in 
km s“^ for [O HI] A5007 and Ha respectively. 

The observed velocity dispersions (uo) - and their Icr 
uncertainties - have been derived from the FWHM measure¬ 
ments of the [O HI]A5007A and Ha lines (cto = 0.4247 x 
FWHM). Corrections for thermal {cth) and instrumental 
{(Ji) broadening have been applied, yielding a final velocity 
dispersion 


a = 





( 1 ) 


with a spin broadening value, in the case of Ha, of 
cr/a(Ha) = 2.4 km s“^ as detailed in Chavez et al. (20141. 


The mean relation between the velocity dispersions of 
[OHI] and of Ha [where we could measure both, including 
our own low-z HHG sample from |Chavez et n[] ( |2014| )] was 
found to be a{Ha) = cr([OHI]) -|- (2.91 ± 0.31) km s“^. We 
used this expresion to estimate cr(Ha) for those objects with 
only cr([OHI]) measured. 


As concluded by|Melnick, Terlevich fc Moles|(|1988|) and 


Chavez et al. (20141, imposing an upper limit to the velocity 
dispersion of logcr(H/?) < 1.8 km s“^, minimizes the prob¬ 
ability of including rotationally supported systems. There¬ 
fore from the observed sample we selected all objects having 
logcr(Ha) < 1.8 km s“^ thus reducing the sample to 6 ob¬ 
jects as indicated in Table column (7). 


The adopted emission line velocity dispersions and their 
1-cr uncertainties are shown in Table column (5). 
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Table 1. Observing log. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

Name 

RA 

s 

date 

slitwidth 

exp. time 


(J2000.0) 


(arcsec) 

VIS, NIR(sec) 

HoyosD2 5 

23 28 41.65 

+00 18 20.0 

29/9/2013 

0.6 

2000, 2400 

Q2343-BM133 

23 46 16.18 

+ 12 48 09.3 

29/9/2013 

0.6 

4000, 4640 

Q2343-BX660 

23 46 29.43 

+ 12 49 45.5 

29/9/2013 

0.6 

8000, 9280 

HoyosD2 1 

23 29 08.20 

+00 20 40.7 

30/9/2013 

0.6 

1600, 2400 

Q2343-BX435 

23 46 26.36 

+ 12 47 55.1 

30/9/2013 

0.6 

4000, 4640 

Q2343-BX418 

23 46 18.57 

+ 12 47 47.4 

30/9/2013 

0.6 

4000, 4640 

Q2343-BX436 

23 46 09.06 

+ 12 47 56.0 

30/9/2013 

0.6 

4000, 4640 

MatsudaC3HAE3 

02 02 37.68 

+01 44 33.2 

30/9/2013 

0.6 

4000, 4640 

HoyosD2 12 

02 28 45.05 

+00 41 32.8 

30/9/2013 

0.6 

3200, 4800 


Table 2. Observed [O III] A5007 and Ha central wavelength and FWHM. 


(1) 

Name 

(2) 

^hel 

(3) 

Ac ([O III] A5007) 
(nm) 

(4) 

FWHM([0 III] A5007) 
(km s“^) 

(5) 

Ac (Ha) 

(nm) 

(6) 

FWHM (Ha) 
(km s“^) 

(7) 

Remarks^ 

HoyosD2 5 

0.6364 

819.359 ± 0.002 

96.6 ± 1.3 

1073.861 ± 0.012 

113.3 ± 7.8 

1 

Q2343-BM133 

1.4774 

1240.440 ± 0.003 

134.1 ± 1.5 

1625.881 ± 0.005 

142.6 ± 2.1 

1 

Q2343-BX660 

2.1735 

1589.181 ± 0.006 

150.9 ± 2.5 

2082.558 ± 0.012 

— 

1 

HoyosD2 1 

0.8510 

926.740 ± 0.015 

118.2 ± 11.4 

— 

— 

1 

Q2343-BX435 

2.1119 

1557.924 ± 0.034 

171.3 ± 15.6 

2042.258 ± 0.021 

178.6 ± 7.2 

2 

Q2343-BX418 

2.3052 

1654.870 ± 0.002 

135.3 ± 1.0 

— 

— 

1 

Q2343-BX436 

2.3277 

1666.438 ± 0.025 

180.0 ± 10.8 

— 

— 

2 

MatsudaC3HAE3 

2.2397 

1622.102 ± 0.054 

292.7 ± 23.4 

— 

— 

2 

HoyosD2 12 

0.6816 

841.988 ± 0.002 

84.4 ± 1.2 

— 

— 

1 


“ 1: Object used in the analysis, 2: log cr > 1.8; object not used. 


3.2 Fluxes 


3.3 Data from the literature 


For the objects selected from Erb et al. ( 2006a|b] ), the Ha 
fluxes were obtained from the literature directly and we 
can readily deduce from the reddening corrected /(Ha), the 
/(H/3) from their theoretical ratio. The three objects from 
[Hoyos et al.| ( |2005| do not have a direct measurement of 
their line fluxes; Their line luminosities were obtained us¬ 
ing a rough estimate from their total blue luminosity and 
equivalent width (see|Terlevich &: Melnick||1981[): 


Hc = -2.51og[L„(H^)/FlTA] + 79.7 , (2) 


where Be is the absolute blue continuum magnitude, 
Lo(H/3) is the reddening corrected H/3 luminosity and EW\ 
is the rest frame equivalent width of H/3. These three ob¬ 
jects have a larger uncertainty in their estimated emission 
line luminosities. 

When available the reddening (here taken as Ay) was 
obtained from the literature. The Av was derived from the 
publisehd E{B — V) using the value of = 4.05 given by 


Calzetti et al. (20001. For those objects where the reddening 


was not available the mean Av = 0.33 from our local sample, 
was adopted. We have verified that the mean values of Av 
for the local and high-z samples are compatible. 

The top block of table column (7) shows the adopted 
H/3 fluxes and their Icr uncertainties, after internal extinc¬ 
tion correction, the adoptedAv values are shown in column 
( 8 ). 


To complement the data and compare results with a larger 
sample (albeit of lower quality velocity dispersions) we have 
searched the literature for measurements of emission line 
FWHM of high-z HHG. Following strictly our selection cri¬ 
teria we were able to select a sample of 6 HHG from |Erb| 
et al. 1 2006a|b' I, 1 from Maseda et al. (20141 and 12 from 


Masters et al. (20141 for which (j(Ha) and /(Ha) are given. 


have line ratios and position in diagnostic diagrams corre¬ 
sponding to HII regions and EW(H/3) > 50A or EW(Ha) > 
200A plus Log (T(Ha) < 1.8 -I- 1 sigma error. Only those 
objects with error less than 25% in the measured velocity 
dispersion are included. 

The lower block of Table |3] shows the data for the 19 
objects selected from the literature. 


4 DISCUSSION 
4.1 The L — a relation 

Taking the concordance ACDM cosmology as the fiducial 
model we calculated the distances and hence the luminosities 
for all the objects in our sample (see subsection 4.3). The 
values of the luminosities are listed in Table column (6). 
As mentioned in section 3.1 three of the objects listed in 
tables 1 and 2 are not included in the analysis because they 
do not fulfil the selection criterion on a. 

Figure [^shows the L — a relation for the 25 high-z sam¬ 
ple of H II galaxies [6 high-z HHG observed with XShooter 










































6 R. Terlevich et al. 


Table 3. Luminosity and gas velocity dispersion of high- 2 : Hii Galaxies obtained from the literature and from observations. 


(1) 


(2) 

(3) 

(4) 


(5) 


(6) 


(7) 



(8) 


(9) 

(10) 

Name 


RA 

5 

^hel 

log 1 

•7 

logL (H/3) 

f (H/1) 




W (Ha) 

Ref. 




(J2000.0) 


(km s 

-u 

(erg s“^) 

(10“11 erg s 

-1 cm-2) 




(A) 


Q2343-BM133 

23 46 

16.18 

+ 12 48 09.31 

1.4774 

1.756 ± 0.017 

42 202"*”^'^^^ 

13.069 ± 

1.293 

0.20 ± 0.14 

2245 

1 

Q2343-BX418 

23 46 

18.57 

+ 12 47 47.38 

2.3052 

1.758 ± 0.016 

42 041'*'^'^^® 

3.031 

± 0.116 

0.14 ± 

0.10 

1639 

1 

Q2343-BX660 

23 46 29.43 

+ 12 49 45.54 

2.1735 

1.808 ± 0.016 

42.024tO;«g 

3.363 ± 0.146 

0.04 ± 0.03 

488 

1 

HoyosD2-5 

23 28 41.65 

+00 18 20.00 

0.6364 

1.597 ± 0.023 

41.393+1^1 

16.171 

± 6.573 

0.28 ± 0.04 

96“ ± 5 

2 

HoyosD2-l 

23 29 08.20 

+00 20 40.70 

0.8510 

1.695 ± 0.049 

41.692+:1®1 

15.803 ± 6.424 

0.28 ± 0.04 

98“ ± 5 

2 

HoyosD2-12 

02 28 45.05 

+00 41 32.80 

0.6816 

1.527 ± 0.027 

41 297'*"^'^^® 

10.965 ± 4.410 

0.00 ± 0.04 

110“ ± 15 

2 

HDF-BX1277 

12 

37 

18.59 

+62 09 55.54 

2.2713 

1.799 

± 

0.062 

4i.yui_Q 056 

2.305 

± 

0.201 

0.29 

± 

0.09 

— 

1 

Q0201-B13 

02 

03 

49.25 

+ 11 36 10.58 

2.1663 

1.792 

± 

0.070 

41.421+;«J9 

0.845 

± 

0.035 

0.01 

± 

0.00 

— 

1 

Q1623-BX215 

16 

25 

33.80 

+26 53 50.66 

2.1814 

1.845 

± 

0.093 

4i.860t“;“i 

2.283 

± 

0.290 

0.41 

± 

0.12 

— 

1 

Q1623-BX453 

16 

25 

50.84 

+26 49 31.40 

2.1816 

1.785 

± 

0.028 


9.076 

± 

2.064 

0.84 

± 

0.25 

187 

1 

Q2346-BX120 

23 

48 

26.30 

+00 20 33.16 

2.2664 

1.792 

± 

0.084 

41 e,+0-042 

41.815_o 046 

1.875 

± 

0.106 

0.02 

± 

0.00 

— 

1 

Q2346-BX405 

23 

48 

21.22 

+00 24 45.46 

2.0300 

1.699 

± 

0.035 

42 125"*”^'^^^ 
44.iZD_o 039 

5.009 

± 

0.083 

0.03 

± 

0.01 

358 

1 

COSMOS-17839 

10 

00 

40.96 

+02 21 38.88 

1.4120 

1.664 

± 

0.084 

41 205"*”^'^^® 

4:i.2UD_j^ g4g 

1.472 

± 

1.446 

0.33 

± 

0.09 

325 ± 230 

3 

WISP159-134 

20 

56 

30.91 

-04 47 56.30 

1.3000 

1.686 

± 

0.045 

41 684'*'^'^^^ 
4H.OO41_0 060 

5.443 

± 

0.532 

0.12 

± 

0.09 

314 ± 36 

4 

WISP173-205 

01 

55 

23.64 

-09 03 10.20 

1.4440 

1.834 

± 

0.045 

4i.684tS:);?^ 

4.196 

± 

0.535 

0.09 

± 

0.15 

603 ± 42 

4 

WISP46-75 

22 

37 

56.48 

-18 42 46.10 

1.5040 

1.839 

± 

0.066 

41.832+;+, 

5.334 

± 

1.794 

0.15 

± 

0.43 

245 ± 28 

4 

WISP22-216 

08 

52 

46.29 

+03 08 45.90 

1.5430 

1.641 

± 

0.040 

4H.ooi_o 062 

3.350 

± 

0.347 

0.00 

± 

0.12 

— 

4 

WISP64-2056 

14 

37 

30.20 

-01 50 51.40 

1.6100 

1.746 

± 

0.039 

41 716+0-050 
^x./xo_0 056 

3.457 

± 

0.304 

0.33 

± 

0.09 

— 

4 

WISP138-173 

15 

45 

31.03 

+09 33 30.00 

2.1580 

1.814 

± 

0.040 

49 11 o+u-Uby 

42.118_o 068 

4.245 

± 

0.504 

0.33 

± 

0.09 

— 

4 

WISP64-210 

14 

37 

28.34 

-01 49 54.40 

2.1770 

1.830 

± 

0.039 

42 043+'^'^^^ 
4IZ.U‘X4_o 059 

3.498 

± 

0.331 

0.33 

± 

0.09 

— 

4 

WISP204-133 

11 

19 

46.37 

+04 10 30.80 

2.1910 

1.765 

± 

0.063 

41 607+’^'^^'^ 

4ix.oui_o 060 

1.262 

± 

0.126 

0.00 

± 

0.00 

— 

4 

WISP70-253 

04 

02 

02.50 

-05 37 19.50 

2.2150 

1.628 

± 

0.041 

41 590"*”*^'*^^® 
d:X.oyu_o 042 

1.182 

± 

0.045 

0.00 

± 

0.00 

— 

4 

WISP96-158 

02 

09 

26.37 

-04 43 29.00 

2.2340 

1.702 

± 

0.043 

41 964+’^''^^^ 
4±x.yO‘X_o 059 

2.742 

± 

0.259 

0.33 

± 

0.09 

— 

4 

WISP138-160 

15 

45 

36.29 

+09 34 26.70 

2.2640 

1.838 

± 

0.044 

42.318+:+, 

5.987 

± 

0.512 

0.33 

± 

0.09 

— 

4 

WISP206-261 

10 

34 

17.56 

-28 30 49.80 

2.3150 

1.693 

± 

0.044 

42-153+; + 

3.880 

± 

1.525 

0.91 

± 

0.49 

— 

4 


NOTE. - Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, 

arcminutes, a nd arcseconds. ____ 

References: (ijErb et aL]j2006b[, (2]Hoyoa et al.|j2005), (SjMaseda et al.||2oT4) & (djlvtastera et al.||2ol4). 


(red stars) and 19 high-z HIIG from the literature (green 
triangles)], and the local sample of GHIIR and HIIG from 
[Ghavez et al. (20141. 

The result is a remarkably tight correlation that 
underpins the use of the L — a relation as a dis¬ 
tance estimator over a wide range of distances, basi¬ 
cally from the local group of galaxies (LMC, SMC, 
NGC 6822, M 33) up to at least z~ 2.3. 

Although here we are only considering the two dimen¬ 
sional L — a relation, we would like to point out that ac¬ 
cording to Chavez et al. (20141 by including additional ob¬ 
servables in the L — a relation like the size of the ionized 
gas region, the equivalent width of either or Ha and the 
ionized gas metallicity or the continuum colour, the scatter 
is substantially reduced from an rms~0.35 to an rms<0.25. 
The importance of this reduction in the scatter of the dis¬ 
tance estimator cannot be overemphasized. 


4.2 The Hubble diagram 

Figure]^ shows the Hubble diagram for the joint sample of 
local and high-z systems. The points correspond to individ¬ 
ual HIIG; their distance moduli are obtained from: 

ti°’’ = 2.5 log L(H^),, - 2.5 log f(H/3) - 100.2 (3) 

where Z/(H/3)o- is estimated from the L—a relation (Figure]^ 
calculated for the joint local HIIG (107 objects) and GEHR 


(24 objects) samples, 

logL(H/3) = (5.05 ±0.097) log a(Ha)-1(33.11 ±0.145) (4) 

The continuous lines show the behviour of the theo¬ 
retical distance modulus with redshift computed for three 
different cosmological models as: 

/r*'* = 51ogDi(p, 2 i) ± 25, (5) 

where Zi) is the luminosity distance as calculated from 

a set of cosmological parameters, p, and a given value of red- 
shift, z. The red line shows the behaviour of the Concordance 
ACDM cosmology with Ho=74.3. The green line shows the 
trend for a cosmology with Q.m = 0.3 and wo = —2.0. The 
solid blue line corresponds to Q,m = 1.0 and Ha = 0.0. In 
all three cases Hfc = 0 i.e Universe flatness is preserved. The 
differential version of the Hubble diagram is shown in the 
bottom panel. 

This is a remarkable and unique Hubble diagram in the 
sense that it covers a huge dynamical range with a single 
distance estimator. It connects galaxies in the local group 
to galaxies at z~ 2.3, a range of almost 30 magnitudes in 
distance modulus or more than 5 dex in redshift. 


4.3 Towards precision cosmology with HII 
galaxies 

The most general set of cosmological parameters, assum¬ 
ing a flat Universe, a negligible value of the radiation den¬ 
sity parameter Hr and using the Chevallier-Polarski-Linder 
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Figure 2 . L — a relation for the combined local (131 HUG and GHIIR) and high-z (25 HUG) samples, the fit corresponds only to the 
local sample of 131 objects. Blue squares: GEHR. Blue dots: local HUG. Red stars: our high-z XShooter observations. Green triangles: 
data from the literature. The inset shows the distribution of the residuals of the fit. The parameters of the fit are indicated at the top. 


(CPL) | |Chevallier fc Polarski|[?001[ |Linder|2003[ ) model for 
parametrizing the value of the dark energy equation of state 
parameter w{z), is given by p = {Jfo, wo, trii} where 
Wo, Wi} are the first two terms of a Taylor expansion around 
the present epoch, namely w{a) = wo + wi(l — a) where 
a = l/(l-|- 2 ) is the scale factor of the Universe. In this 
case the luminosity distance 
(P : Zi), is given by (cf. e.g. 

2003]), 


used to estimate the value of 


Weinberg 2008 Frieman et al. 


Dl 



Hq + z)^ + (1 


X exp 


/ —Swiz\ 

U+i / 


( 6 ) 




(7) 


To restrict the set of cosmological parameters we min¬ 


imised the Likelihood function, 

2/ >, fi) — 

X (pj - -“2-> (°) 

i=l 

fi) are the ‘observed’ distance moduli obtained from 
equation!^ ai are the measured velocity dispersions and fi 
are the measured H/3 fluxes for each object. Zi) are the 

‘theoretical’ distance moduli from equation [^obtained from 
the measured redshifts by using a particular set of cosmo¬ 
logical parameters p. o^ob are their errors propagated from 
the uncertainties in at and fi and the slope and intercept 
of the distance estimator in equation The summation is 
over the combined sample of HIIGs, 

Adopting the value of Ho = 74.3 ± 3.1 obtained in 
[Chavez et al.| ( |2012| ), Wi = 0 we obtain the results shown 
in Figure [^ for p = {flm,wo} for the high-z sample only. 
Panel (a) shows the solution for the 6 high-z HIIG observed 
with XShooter; Panel (b) shows the effect of including the 
19 high-z objects from the literature. It is encouraging that 
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10"^ lO"* 10"^ 10"^ 10"' 10° 

z 

Figure 3. Hubble diagram for our sample of low and high-z HUG for three different cosmologies. The solid red line indicates the 
concordance ACDM cosmology with flrn = 0.3; uiq = —1.0 and Ho=74.3. The solid green line shows a cosmology with Qrn = 0.3 and 
wq = —2.0. The solid blue line corresponds to = 1-0 and Qx = 0.0. In all three cases Qj. = 0. Residuals are plotted in the bottom 
panel. Note the huge dynamical range in distance modulus of almost 30 magnitudes covered with the L — a distance estimator. 


even for this small number of high- 2 ; objects we are able to 
restrict so much the solution space. Although the literature 
data are of lower quality regarding the FWHM measure¬ 
ments, still by adding them we are able to improve the result. 
This is a remarkable result considering the small number of 
intermediate to high -2 data (just 25 objects), of which only 
6 have high quality ad-hoc observations. 

In Figure we compare our results for the space p = 
{nm,rco}, joining the high -2 with the local HIIG samples 
(see left panel), with recent results from SNe la, CMB and 
BAO (right panel). The figure shows the contrains of the 
properties of dark energy using SNe la alone | |Amanullah| 
et al.|20l0 1, the seven-year Wilkinson Microwave Anisotropy 


Probe data of the CMB ( Komatsu et al.|2011 1, the position 
of the BAO peak from the combined analysis of the SDSS 
Data Release 7 and 2dFGRS data ( Percival et al.|2010 1. The 
combined restrictions from SNe la, CMB and BAO and the 
measurement of the Hubble constant (Ho) from Cepheids 
i Riess et al.|20ll| are also shown. 

It is clear from the figure that our present constraints 
on the space p = are weaker than those for SNe 

la, but this is not surprising since in our case we have only 


156 objects, most of them at z<0.16 a region of space where 
differences between Cosmological models are almost negli¬ 
gible, vs. 580 SNe la with a maximum redshift of ~ 1.5. 
The strength of our results is that our sample includes 19 
objects with 2 > 1.5 where the differences between models 
reach maximum values. 

From the comparison of the figures we can conclude first 
that there are no systematic shifts between the HII galaxies 
and SNe la solutions and secondly that with a larger sample 
of HII galaxies with high quality data it may be possible 
to achieve at least similar and probably even better results 
to those obtained with SNe la as found in our simulations 
( [Plionis et al.|[^ll[ ) and discussed below. 


4.4 The future of the Z/(H/3) — a distance estimator 


To estimate the number of high -2 tracers (Nhz) required to 
substantially reduce the cosmological parameters solution 


space we used the figure of merit (FoM) as in Plionis et al. 
( |2011[ ). The FoM is the reciprocal area of the 2cr contour 
in the parameter space of any two degenerate cosmologi- 
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Om 

(a) (b) 

Figure 4. Solution space in the plane {r2,„,uio} (see text). Panel (a) for the 6 XShooter high-z objects. Panel (b) Same as panel (a) 
including 19 high-z objects from the literature. In both panels we show the 1 and 2 (t contours. 



Figure 5. Comparison of restrictions on the plane {Qm, rco}(see text). Panel (a) shows our results, obtained as described in the text 
for the combined 25 high-z HUG and the local sample (131 HUG and GHIIR). 1 and 2cr contours (random) are shown. Panel (b) after 
[Suzuki et al.|)[2012|l shows the recent results for 580 SNe la, CMB and BAOs, the 1, 2 and 3cr contours (random) are shown. 


cal parameters. In this way a larger FoM indicates better 
restrictions to the cosmological parameters. jPlionis et al.j 
(20111 use the parameter S or “reduction factor” to com¬ 
pare the ratio of FoM of SNe la -I- high-z tracers to that of 
only SNe la. They found that the number of high-z tracers 


needed to obtain a given S in the quintessence dark energy 
(QDE) model implying that w is constant but different from 
-1, when combined with the intermediate and low-z SNe la, 
can be expressed as: 


Nh. ~ ISTS/S-ioo - 88 


(9) 
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where Sioo = 1.871og((cr^)~^ +0.74) + 1.28, and (cr^) is the 
mean distance modulus error for the tracer. In the case of 
HUG, (T^ ~ 0.6, 5100 — 1.99 and consequently in order to 
obtain a reduction of 5 = 2, we need around 100 high-z 
HUG. We can visualize this in Figure panel (a) for the 
p = {f2m,rco} plane, where a factor of 2 reduction would 
mean that the future 2(j contours would be similar to the 
present Icr ones. 

For the GPL model where wi is variable the same result 
can be expressed as: 

Nh, ~ 4045/Sioo - 300 (10) 

where 5ioo = 0.491og((cr^)~^ + 0.65) + 1.09. In the case 
of HUG 5ioo — 1.27 and consequently in order to obtain a 
reduction of 5 = 2, about 300 high -2 HIIG would be required 
to achieve a result similar to that obtained with present day 
SN la samples. 


5 CONCLUSIONS 

We have tested the use of high redshift HIIG to trace the 
expansion of the Universe by means of their I/(H/3) — a dis¬ 
tance estimator. To this end, we presented observations of a 
sample of only 9 HIIG in the redshift range of0.65S2^2.3 
obtained with the ESO VLT XShooter spectrograph. After 
rejecting three HIIG due to either poor data or because they 
fall outside the selection window, we have used the remain¬ 
ing six to obtain constraints on the {ilo, Hm} and {flm, uio} 
planes. 

The results are surprisingly good considering the small 
number of objects and that only two nights of XShooter time 
were used. The addition of 19 HIIG from the literature (al¬ 
beit of poorer quality) bringing the total number of objects 
to 25, did improve the results and provides a clear indica¬ 
tion of how much an increase in number can improve the 
constraints on the {77o,fIm}, {IIm,wo} and{wo,Wi} planes. 

Although our constraints are consistent with other de¬ 
terminations they are, as expected, definitely weaker. This is 
due to the small size of the intermediate to high -2 sample, 
and the considerable uncertainties in the data taken from 
the literature. 

Using the figure of merit approach we have estimated 
the expected improvement in the estimates of the cosmo¬ 
logical parameters with a larger sample of high -2 HIIG. In 
particular, between 100 and 300 high -2 HIIG are needed to 
obtain a reduction in the errors of at least a factor of two. A 
factor of 2 reduction implies that the 2cr contours in Figurej^ 
panel (a) for the p = {HmjWo} plane would be similar to 
the present la ones and comparable to the 2a SNe la con¬ 
tours in panel (b) that are the result of the analysis of 580 
SNe la that took at least a decade and hundreds of nights of 
large telescope time to compile. In contrast, the observation 
of 100 HIIG can easily be achieved with multiple integrated 
field units instruments like VLT- KMOS, in less than thirty 
hours observing given the relative abundance of HIIG with 
a considerable number of them appearing in its 7.5 arcmin 
field of view. 

Finally the comparison of the results of SNe la and HIIG 
will undoubtedly contribute to learning about the system¬ 
atic errors that limit the precision of both empirical methods 


while helping us to gain an insight into the intrinsic prop¬ 
erties of HIIG at high redshifts. We also envisage that a 
substantial improvement to the present restrictions of cos¬ 
mological parameters will be obtained by combining a few 
hundred HIIG with the SNe la data. 
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